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Abstract. 

An exact quantization rule for the Schrodinger equation is presented. In the exact quanti- 
zation rule, in addition to Nit, there is an integral term, called the quantum correction. For 
the exactly solvable systems we find that the quantum correction is an invariant, independent 
of the number of nodes in the wave function. In those systems, the energy levels of all the 
bound states can be easily calculated from the exact quantization rule and the solution for the 
ground state, which can be obtained by solving the Riccati equation. With this new method, 
we re-calculate the energy levels for the one-dimensional systems with a finite square well, with 
the Morse potential, with the symmetric and asymmetric Rosen-Morse potentials, and with the 
first and the second Poschl- Teller potentials, for the harmonic oscillators both in one dimension 
and in three dimensions, and for the hydrogen atom. 



In the development of quantum mechanics, the Bohr-Sommerfeld quantization rules of 
the old quantum theory [1] occupy a position intermediate between classical and quantum 
mechanics. The WKB approximation [2-4] is a method for the approximate treatment of 
the Schrodinger wave function with another quantization rule [1]. Both quantization rules 
are approximate. We are going to derive an exact quantization rule for the one-dimensional 
Schrodinger equation: 

^ {x) = "IP 

where M is the mass of the particle, and the potential V(x) is a piecewise continuous real 
function of x satisfying for definiteness 



—^(x) = -—[E-V(x)]^(x), (1) 



V(x) > E, — oo < x < xa or xb < x < oo, 

V(x) = E, x = xa or x = xb, (2) 

V(x) < E, xa < x < xb, 

where xa and xb are two turning points. Between two turning points, the momentum is 
k(x) = y/2M[E-V(x)]/h. 
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Yang pointed out in a talk on monopole: "For the Sturm-Liouville problem, the fundamen- 
tal trick is the definition of a phase angle which is monotonic with respect to the energy." [5] 
This phase angle is the logarithmic derivative <j>(x) = tp(x)~ 1 dtp(x)/dx of the wave function 
tp(x). Due to the Sturm-Liouville theorem, <f>(x) at any given point x = xq is monotonic with 
respect to the energy. From the Schrodingcr equation (1), <f>(x) satisfies the Riccati equation: 



d A( \ 2M [ P 

Tx m = -1F [E 



V(x)]-<j>(x) 2 . 



(3) 



It shows that 4>{x) decreases monotonically with respect to x between two turning points 
where E > V(x). Note that as x increases across a node of the wave function ip{x) where 
E > V(x), 4>{x) decreases to — oo, jumps to +oo, and then, decreases again. The Riccati 
equation is a differential equation of the first order so that it is much easier to find a special 
solution from the Riccati equation than from the Schrodinger equation. 

Letting tan#(x) = k(x)/(p(x), we have 6{x) = Arctan [k(x)/cp(x)] +nir, where Arctan (3 
denotes the principle value of the inverse tangent function, — ir/2 < Arctan j3 < tt/2, and n 
increases by one as x increases across a node of <f>(x) where E > V(x). Then, we have 



d0(x) 



( Hx) \ 



dx = Ntt — lim Arctan [ 7 I + lim 

dx x^x A + \(f>(XAj J x^x B - 



Arctan 



k{x) 
<j>(x B ) 



(4) 



where N is the number of nodes of 4>{x), and two terms with limit arc vanishing if the potential 
V(x) is continuous at the turning points. From Eq. (3) we have 



d9(x) 
dx 



k(x) — <p(x) 



dk(x) 
dx 



d(f)(x) 



-i -l 



dx 



E > V(x). 



(5) 



Integrating both sides of Eq. (5) over the interval between two turning points, we obtain the 
quantization rule without any approximation: 



f 



k(x)dx — Ntt - 



f 



<f>(x) 



dk(x) 



dx 



d<p(x) 



dx 



dx. 



(6) 



Since 4>{xa) > 0, 4>{xb) < 0, and <j>(x) decreases monotonically in the region xa < x < xb, 
the number N of nodes of 4>(x) in that region is larger by one than the number of nodes of 
the wave solution ip(x). Due to the Sturm theorem, the number (N — 1) of nodes of ip(x) 
increases as the energy E increases. 

Generalize the quantization rule to the three-dimensional Schrodinger equation with a 
spherically symmetric potential. After separation of the angular part of the wave function, 
i)j(r) = r~ 1 R(r)Ym(0, <p), the radial equation of the Schrodinger equation is 



d 2 R(r) 2M r , „ „, , 



U{r) = 



[£ + i)h 2 

2Mr 2 



+ V(r). 



(7) 



Since Eq. (7) is similar to Eq. (1), the quantization rule (6) is generalized to the three- 
dimensional Schrodingcr equation with a spherically symmetric potential 



rrB rrB 

/ k(r)dr = Ntt + <f>(r) 

J v a J r A 



dk(r) 
dr 



d<p{r) 
dr 



dr. 



(8) 



The quantization rule (6) [or (8)] is proved without any approximation, so that it is exact. 
The first term Ntt in the quantization rule relates to the contribution from the nodes of 
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the wave function, and the second term is called the quantum correction. We find that 
the quantum correction is independent of the number of nodes of the wave function for the 
exactly solvable systems. In those systems, one is able to calculate easily the energy levels of 
all bound states from the exact quantization rule and the solution of the ground state, which 
can be calculated directly from the Riccati equation. With this new method, we are going to 
re-calculate the energy levels of those systems. 

The one-dimensional system with a finite square well, where V(x) = V A when x < —n, 
V(x) — Vb when x > ir, and V(x) = when — ir < x < n, is a typical example in the course of 
quantum mechanics. However, the solutions can be obtained with the exact quantization rule 
much simpler than with the standard method. Since k(x) is constant between two turning 
points, the quantum correction in the quantization rule is vanishing. For the finite square 
well, two terms with limit in Eq. (4) has to be included in the quantization rule (6) because 
the potential jumps at the turning points. The energy levels of the system is directly given in 
the quantization rule (6) and (4) [6]: 



2nk N = Ntt — Arctan ( ' /V ) + Arctan ( — ) , (9) 



(xa)J \4>n(xb) 



k N 



where (f> N (x A ) = ^2M(V A - E N )/h, 4>n{x b ) = -y/2M(V B - E N )/h. When V A = V B ^ oo, 
we have k N = N/2 [1]. 

The potential for a one-dimensional harmonic oscillator is V(x) = Mui 2 x 2 /2. The turning 
points are solved to be x B = —x A = a^ 1 ^2E n /(hui), where a — y/Mcj/h, and (n + 1) 
is the number of nodes of <f) n {x). The momentum between two turning points is k n (x) = 
a 2 [(x — x A ) (x B — x)] 1//2 . The solution with one node and no pole only has the form as 
4>o(x) = —Cx where C > due to the monotonic property. Substituting it into the Riccati 
equation (3) with the potential of a harmonic oscillator, we find <fio(x) = —a 2 x with Eq = 
hu/2. Evidently, <Pq{x) is negative when x — > oo and is positive when x — ► — oo, so that 
the solution satisfies the physically admissible boundary conditions. All the solutions in the 
following examples have the similar behavior. Two integrals in the exact quantization rule (6) 
are calculated to be 



p \dk (x) 

Jxa 



L 



dx 

x B 



d4> (xy 



dx = -tt/2. (10) 



dx 

k n {x)dx = E n ir/(huj). (11) 

'X A 

The quantization rule (6) coincides with the quantization rule in WKB approximation: 

k n (x)dx = (n + 1/2)tt. (12) 



/ 

'J X t 



The energy levels for the one-dimensional harmonic oscillator are [1]: 

E n = {n + l/2)huj. (13) 

The one-dimensional Morse potential is V(x) — D ( e _2;r / a — 2e~ x / a ) = Dy(y — 2), where 
y = e~ x / a . The turning points are y A = e~ XA l a = 1 — yjl + E n /D and y B = e~ XB l a = 
l + yT + E n /D, where (n+1) denotes the number of nodes of the logarithmic derivative <j> n (x). 

1/2 

The momentum between two turning points is k n (x) = [(2MD/h 2 )(y — y A )(y B — y)] ■ The 
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solution with only one node has to be <j>o(x) — C1J/ + C2 with C\ > 0. Substituting it into the 
Riccati equation (3) with the Morse potential, we find <f>o(x) = V2MD(y-l)/h+l/(2a) with 
E a = — [y/~D — h/ (a\/8M)] 2 . Two integrals in the exact quantization rule (6) are calculated 
to be 



<j>o(x) 



dk (x) 



dx 



d<j> (x) 



dx 



dx 



-tt/2, 



k n (x)dx 



XA 

The quantization rule is f*^ k n {x)dx = (n + l/2)n. Thus, the energy levels for the one- 
dimensional system with the Morse potential are [7] : 



a\/2M 



VD- ^T-E~ n 



(14) 



(15) 



E„ — — 



Vd- 



(2n + l)h 
2aV2M 



(16) 



The asymmetric Rosen-Morse potential [7] in one dimension is V(x) — ~Uosech 2 (x/a) + 
Ui tanh(x/a), where < Ui < 2Uq. If U\ — 0, V(x) is called the symmetric Rosen-Morse 
potential. Let y = tanh(x/o), yA — tanh(x a I a ) , an d Vb = tanh(a; B /a), where xa and xb are 
two turning points satisfying V(xa) — V(xb) = E n . We have 



VA = 



El 

2U () 



2U 



E„ 

— + 1, 

U 



VB = 



El 

2U () 



El 

2U Q 



U + ' 



where (n + 1) is the number of nodes of the logarithm derivative (f> n (x). The momentum 

between two turning points is k n (x) = \J 2MUo/h 2 [(y — jmXj/b — y)} 1 ^ 2 - From the Riccati 
equation (3) with the asymmetric Rosen-Morse potential we obtain the solution with one node 



Mx) = ~2a 



1 



8a 2 MU 

h 2 



1/2 



-1 }y 



MUi 



2 n 2 



with E 



MU 2 



2M 2h 2 C 2 



Now, two integrals in the quantization rule (6) are calculated to be 

1 n -1 



rxi 
Jx A 

f 



<f>o{x) 



dk (x) 



dx 



d(j> (x) 



dx 



dx 



k n (x)dx 



air y/2MU 



2 n 

The quantization rule (6) becomes: 

1 



2- 



-E n - U x 



nc 



U 



-En + Ul 

U 



+ 



C/ ' 2 
Thus, the energy E n is [7] 



E n + Ux _ {n+l)h y/2MU _ (Ca - n)% 

ay/2MU 



U 



+ 



ay/2MU nC 



Ti 2 {C-n/af 
2M 



MU 2 



2h z {C-n/a) 2 



(17) 
(18) 

(19) 

(20) 
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The condition of existence for the bound state whose wave function has n nodes is U\ < 
h 2 (C — n/a) 2 /M < 2Uq. When U\ = 0, the asymmetric Rosen-Morse potential becomes the 
symmetric one. The energy levels (20) with U\ = hold for the symmetric Rosen-Morse 
potential [7]. 

The first Poschl- Teller potential [7] in one dimension is 



V(x) 



n 2 



2Ma 2 



M/x-l) , A(A-l) 



sin (x/a) cos 2 (x/a) 



an 

0<x<—, 



(21) 



where n and A are constant greater than one. The potential V(x) tends to infinity as x tends 
to or air/2. Let y = tan 2 (x/a), ua = t&n 2 (xA/a), and ys = tan 2 (a;_e/a), where xa and xb 
are two turning points satisfying V(xa) — V(xb) = E n . We have 



VA + VB 



2Ma 2 E n /h 2 - fi(fi - 1) - A(A - 1) 

x(x^T) ■ 



VAVB 



A(A-l)' 



where (n+ 1) denotes the number of nodes of the logarithm derivative <j) n {x). The momentum 
k n (x) between two turning points is k n {x) = [A(A~ l)(y — VA){yB — J/)/j/] 1 ^ 2 /a when E > V(x). 
The solution of the Riccati equation with only one node is (f>o(x) = —\y 1 / 2 /a + [iy~ 1 / 2 /a with 
E = h 2 (p + X) 2 /(2Ma 2 ). 

Two integrals in the quantization rule (6) are calculated to be 



0o (x) 



dko(x) 



dx 



d(j> (x) 



dx 



dx = — 
2 



(p + X-2)- y/fl(fj, - 1) - y/\(\ - 1) 



k n (x)dx = — 



ay/2ME n 



y/ -l)-y/\(\-l) 



(22) 



(23) 



The quantization rule (6) becomes ay / 2ME n /h = 2(n + l) + (fi + A — 2), from which we obtain 
the energy E n [7] 

h 2 ( f i + \ + 2n) 2 



E„ 



2Ma 2 



The second Poschl- Teller potential [7] in one dimension is 



V(x) 



2Ma 2 



Mm-i) a(a + i) 



sinh (x/a) cosh (x/a) 



(24) 



(25) 



where A > fi — 1 > 0. The potential V(x) tends to infinity at x = 0. Let y — tanh 2 (x/a), 
yA = tanh 2 (xa /a) , and ys — tanh 2 (xs/a) where xa and xb are two turning points where 
V(x A ) = V(x B ) = E n . We have 



VA + VB 



2Ma 2 E n /h 2 + nifi - 1) + A(A + 1) 
A(A + 1) : 



VaUb 



A(A + 1)' 



where (n+ 1) denotes the number of nodes of the logarithm derivative <p n (x). The momentum 
k n (x) between two turning points is k n (x) = [A(A+l)(y— yA)(yB— y)/!/} 1 ^ 2 /&, when E > V(x). 
The solution of the Riccati equation with only one node is (f>o(x) = —\y 1 / 2 /a + [iy~ 1 / 2 /a with 
E = -h 2 (X- fi) 2 /(2Ma 2 ). 

Two integrals in the quantization rule (6) are calculated to be 



(j) (x) 



dko(x) 
dx 



d(f>o(x) 
dx 



dx = 



(/l - A - 2) - y/(i(fl - 1) + y/X(X + 1) 



(26) 
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/' 



k n (x)dx = 



a^-2ME n 
h 



- y/pfa - 1) + y/X(X + 1) 



(27) 



The quantization rule (6) reads a^/~2ME n /h — 2(n + 1) + {fi — X — 2), and the energy E n 
is [7] 



En — — 



n\\-n-2nf 



0<n< (X-fj,)/2. 



(28) 



2Ma 2 

The effective potential for the three-dimensional harmonic oscillator is Ui(r) = £(£ + 
l)h 2 /(2Mr 2 ) + Muj 2 r 2 /2. The turning points arc r A = a' 1 {[E ne / {Tiuj)] - [(E ni /(?uj)) 2 - 
£( £ + l)] 1 / 2 } 1 / 2 and r B = ^{[^/(M] + [(£W(M) 2 - ^ + 1)] 1/2 } 1/2 , where a = 
y Mu>/h, (n — £) is a non-negative even integer, and (n — ^ + 2)/2 is the number of nodes 
of the logarithmic derivative (j) n i{r). The momentum between two turning points is k n t(r) = 
(a 2 /r)[(r 2 — r 2 A ){r 2 B -r 2 )] 1 / 2 . The solution with only one node has to be 4> n e(r) = C\r~ Y +C2r 
where C\ > and Ci < 0. Substituting it into the Riccati equation (3) with the harmonic 
oscillator potential, we find <f> n i(r) = {£+ l)r _1 — a 2 r with E n i = hui(£ + 3/2) and n = £. 

Two integrals in the quantization rule (8) are calculated to be 



•Ir a 



<p n e(r) (dk ne (r)/dr) 



Ki{r) 



dr 



dr 



£- y/£(£+ 1) - 1/2 tt/2, n = £. (29) 



/' 

Jta 



k ni (r)dr = 



VW+T) 



The quantization rule (8) becomes J^ B k n i{r)dr 



tt/2. (30) 
T) + 3/2]7r/2. Thus, the energy 



levels for the three-dimensional harmonic oscillator are [1]: 

E n e = (n + 3/2)huj. 

The effective potential for the hydrogen atom is Ug(r) — £(£ 
£ > 0, the turning points r A and tb satisfying UeX r A) = U^(tb) 



l)Ti 2 /{2Mr 2 ) 
= E ne are 



r A = (-2E nl )- 1 [e 2 - [e 4 + 
r B = {-2E n ?y 1 [e 2 +[e i + 



2£(£+l)h 2 E ne /M] 1/2 }, 
2£(£+l)h 2 E nt /M] 1/2 }, 



(31) 

e 2 /r. When 



(32) 



where we denote by (n — £) the number of nodes of the logarithmic derivative (f> n e{r)- When 
£ = 0, we define ta — with U(r A ) = —e 2 /rA ~ — oo, and tb is solved from U(tb) = E n i- 
Thus, equation (32) still holds for £ = 0. The momentum between two turning points is 

k n e(r) = (ft?") -1 [-2Mi? n f(r — ta){tb — f)] 1 ^ 2 ■ Fot £ = 0, near the origin wc have k n g(r) <~ 
r -1 / 2 , and 4> n i(r) ~ r -1 , so that the limit terms in Eq. (4) still vanish. The solution with 
only one node has the form as 4> n t{r) = C\r~ x + Ci where C\ > 0. Combining it with the 
Riccati equation (3) for the hydrogen atom, we find (f> n e(r) = (£ + l)/r — Me 2 /[{£+ l)h 2 ] with 
E ne = -Me A /[2(£ + Ifh 2 ] and n = £ + 1. 

Two integrals in the quantization rule (8) are calculated to be 



Jta 



<t>nt{r) 



dk„i{r) 



dr 



Kt(r) 



dr 



dr 



£ - ^J£{£ + 1) 



£+1. 



f 



k nt (r)dr = 



M 

-2E ne 



-y/£(£+l) 



71". 



(33) 



(34) 
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The quantization rule (8) becomes k n t(r)dr — [n — \J i{l + l)]7r. Thus, the energy levels 
for the hydrogen atom are [1]: 

Me 4 

In this Letter we present an exact quantization rule (6) for the one-dimensional Schrodinger 
equation and (8) for the three-dimensional Schrodinger equation with a spherically symmetric 
potential. We find that the quantum correction term is independent of the number of nodes in 
the wave function for the exactly solvable quantum systems. In such cases, the energy levels 
of the quantum system can be easily solved from the exact quantization rule and the solution 
of the ground state calculated directly from the Riccati equation. For the non-exactly solvable 
systems, one can use the series form of the quantization rule for numerical calculation [8]. 

As is well known, the wave functions and the energy levels for the exactly solvable systems 
can also be solved by the supersymmetry in the recursive way [9] . The logarithmic derivatives 
4>{x) are nothing but the superpotentials in the supersymmetric quantum mechanics. The 
shape invariance of the superpotentials seems to be related to the condition of the quantum 
correction being an invariant. As far as the wave function ip^{x) with (N — 1) nodes is 
concerned, we prefer to calculate its logarithmic derivatives </>n(x) directly from the Riccati 
equation (3). Although there are more or less differences for different examples, the solution 
for 4>n(x) can be taken as a fraction, where the numerator is a polynomial of order N and the 
denominator is a polynomial of order (N — 1), because 4>n(x) contains N nodes and (N — 1) 
poles. Without loss of generality, the coefficient of x N ~ x in the denominator can be chosen 
to be one. Substituting it into the Riccati equation (3), one obtains a coupled algebraic 
equations of order two for the coefficients. Therefore, 4>n(x) can be solved. We have solved 
some examples for the lower excited states. 
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